from sympy import *

x=Symbol('x')
ans=poly((1+x**2+x**4+x+x**3+x**5)*(1+x**5+x**10+x**15))
print(ans)